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Abstract 

We study the persistent current and the Drude weight of a system 
of spinless fermions, with repulsive interactions and a hopping impurity, 
on a mesoscopic ring pierced by a magnetic flux, using a Density Matrix 
Renormalization Group algorithm for complex fields. Both the Luttinger 
Liquid (LL) and the Charge Density Wave (CDW) phases of the system 
are considered. Under a Jordan- Wigner transformation, the system is 
equivalent to a spin- 1/2 XXZ chain with a weakened exchange coupling. 
We find that the persistent current changes from an algebraic to an ex- 
ponential decay with the system size, as the system crosses from the LL 
to the CDW phase with increasing interaction U. We also find that in 
the interacting system the persistent current is invariant under the impu- 
rity transformation p — * 1/p, for large system sizes, where p is the defect 
strength. The persistent current exhibits a decay that is in agreement 
with the behavior obtained for the Drude weight. We find that in the LL 
phase the Drude weight decreases algebraically with the number of lattice 
sites N, due to the interplay of the electron interaction with the impurity, 
while in the CDW phase it decreases exponentially, defining a localization 
length which decreases with increasing interaction and impurity strength. 
Our results show that the impurity and the interactions always decrease 
the persistent current, and imply that the Drude weight vanishes in the 
limit N — > oo, in both phases. 
PACS: 71.10.Pm, 73.23.Ra, 73.63.-b 
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1 INTRODUCTION 



The experimental discovery of persistent currents in mesoscopic rings pierced by 
a magnetic flux, 1-3 earlier proposed theoretically, has revealed interesting new 
effects. The currents measured in metallic and semiconducting rings, cither in 
a single ring or an array of many rings, generally exhibit an unexpectedly large 
amplitude, i.e., larger by at least one order of magnitude, than predicted by 
theoretical studies of electron models with either disorder or electron-electron 
interaction treated perturbatively. 5,6 It has been suggested that the interac- 
tions and their interplay with disorder are possibly responsible for the large 
currents observed, expecting that the effect of the interactions could counter- 
act the disorder effect. However, no consensus has yet been reached on the 
role of the interactions. In order to gain theoretical insight, it is desirable to 
perform numerical calculations which allow to consider both interactions and 
disorder directly in systems with sizes varying from small to large. Analytical 
calculations usually involve approximations which mainly provide the leading 
behavior of the properties for large system sizes. Persistent currents in meso- 
scopic rings strongly depend on the system size, since they emerge from the 
coherence of the electrons across the entire system. Hence, it is most important 
to study the size dependence of the current beyond leading order in microscopic 
models, for a complete understanding of the experimental results. Exact diag- 
onalization was used to calculate persistent currents in systems with very few 
lattice sites. 7 ' 8 In this work, we use the Density Matrix Renormalization Group 
(DMRG) algorithm, 9-11 to study a simplified model incorporating interactions 
and a single impurity, accounting for disorder, in larger system sizes. We con- 
sider a system of interacting spinless electrons on a one-dimensional ring, with a 
single impurity, and penetrated by a magnetic field. We study an intermediate 
range of system sizes, where analytical results obtained by bosonization tech- 
niques for large system sizes, do not yet fully apply. Without impurity, and at 
half-filling, the system undergoes a metal-insulator transition from a Luttingcr 
Liquid (LL) 12 to a Charge Density Wave (CDW) 13 groundstate. The persistent 
current of the interacting system with an impurity was studied before with the 
DMRG, in the LL phase. 14 Here we study the persistent current, and also the 
Drudc weight characterizing the conducting properties of the system, in both 
the LL and the CDW phase, investigating the interplay between the impurity 
and the interactions in the two phases. In mesoscopic systems the separation 
between metallic and insulating behavior is not always obvious, since the local- 
ization length can be of the order or significantly larger than the system size. 
Hence, a finite Drude weight and a current can be observed in the CDW phase 
of a mesoscopic system. It is therefore of great interest to characterize the per- 
sistent current and the Drude weight in both the LL and the CDW phases of 
mesoscopic systems. Although the simple model that we consider is not the 
most appropriate to describe the experimental situation, we hope to obtain use- 
ful information for the understanding of the more realistic systems. Under a 
Jordan- Wigner transformation, 15 the system considered is equivalent to a spin- 
1/2 XXZ chain with a weakened exchange coupling. Hence, our results also 
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provide insight into the spin transport in this type of systems. 



2 THE MODEL 

The Hamiltonian describing a system of spinlcss fcrmions on a ring pierced by 
a magnetic flux, with repulsive interactions and a single hopping impurity, or 
defect, is given by, 

H = H t +Hu, (1) 

where 

N 

H t = -tJ2 (e lHN c]c 1+1 + e -^ C t +1 c,)+(l-p)i (e^cU + e^ N c\c N ) 

3 = 1 

(2) 

is the hopping term, <j> = 27r$/$ Q contains the magnetic flux <3? in units of the 
flux quantum <I> = hc/e, p measures the strength of the defect with values 
between and 1, (p = 1 corresponding to the defectless case), and 

N 

H U = U ^2 n 3 n 3+l ( 3 ) 

3 = 1 

is the interaction term, with U > representing the nearest neighbor Coulomb 
repulsion, and rij = c^Cj, where cj and Cj arc the spinless fermion operators 
acting on the site j of the ring. We consider a system of N sites, with TV even, 
and at half-filling, when M = N/2 particles are present. The lattice constant is 
set to one and periodic boundary conditions, c/v + i = ci, are used. 

Via the gauge transformation Cj — > e~ l ^/ N Cj, the flux can be removed from 
the Hamiltonian, but in the impurity term where the flux is trapped, and the 
quantum phase <j> is encoded in a twisted boundary condition cn+i = eT l ^c\. It 
is then clear that the energy is periodic in (j> with period 2ir, i.e., it is periodic 
in the flux <!> threading the ring with period $o- 16 

After a Jordan- Wigner transformation, Eqs. (2) and (3) can be rewritten, 
respectively, as 

7 N T 
H J = "a E ( S 3 +S 7+i + S I S t+i) + C 1 - P) 2 ^ S t S N + ^SiSn) (4) 

3=1 



and 



N 

H A = AJ2S*S* +1 (5) 

3=1 

with t = J/2 and U = A, and the boundary conditions S^ +1 = (-l) M+1 e^5 ] f 
and Sjv + i = . Hence, the model (1) of spinless fermions is equivalent to a 
spin-1/2 XXZ chain with a weakened exchange coupling, and twisted boundary 
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conditions in the transverse direction. The half-filled case corresponds to total 
spin projection S z = 0. 

The persistent current generated on a ring pierced by a magnetic flux, at 
temperature T = 0, can be obtained from the ground state energy E (<f>), by 
taking the derivative with respect to 0, 

m = -^. (6) 

For the spinless fermion system, Eqs. (2) and (3), I(4>) corresponds to the 
ground state value of the charge current operator I c = it X^jLi ( C ] C J +1 — c j+i c i) ' 
while for the XXZ chain, Eqs. (4) and (5), it corresponds to the ground state 
value of the spin current operator f s = ijJ2f=i {Sj^Sj+i — Sj~ +1 Sy) ■ As a 
consequence of the periodicity of the energy, the current is also periodic in 0, 
with period 2n. Hence, it can be expressed as a Fourier series, 

oo 

J(0)=£j fc sin(fc$. (7) 
fc=i 

and the behavior of the current can be analyzed in terms of the coefficients Ik- 14 
In the noninteracting case (U = 0), it has been found that for large system 
sizes, the current is invariant under the defect transformation p — ► 1/p, 17 ' 18 i.e., 

I(p)=I(l/p). (8) 

We shall investigate the existence of this kind of invariance in the interacting 
case (U > 0), both in the LL and the CDW phases. 

The Drude weight was proposed by Kohn as a relevant quantity to distin- 
guish between a metal and an insulator. 19 It is defined as 

D = N^\, = , m , (9) 

where 4> m is the location of the minimum of E {<j>), which depends on the parity 
of the number of electrons, i.e., (f>. m = or ir for, respectively, an odd or an even 
number of electrons. For the spinless fermion system D represents the charge- 
stiffness and measures the inverse of the effective mass of the charge carriers. 20 
In a metallic conductor D tends to a finite value whereas in an insulator D 
vanishes with the system size N, when N — > oo. In the insulating state the 
Drude weight decays as D ~ exp(— N/£), where £ measures the localization 
length. For the XXZ chain the Drude weight represents the spin-stiffness. 

In a model of free fermions (U — 0) with no impurity (p = 1), it is straight- 
forward to see that the leading behavior of the persistent current I((j>) in the 
system size N, has a saw-tooth like shape with slope —vf/(itN), where vp is the 
Fermi velocity. Thus, the amplitude of the current scales with 1/N, vanishing 
in the limit TV — > oo. The discontinuity in /(</>), that results from a degeneracy 
of energy levels associated to the translation symmetry, 8 appears at <fi = or 
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7r for, respectively, an even or an odd number of electrons. In the presence of 
an impurity (p ^ 1), bosonization 21 and conformal field theory 17 calculations 
predict that the shape of the current I{<j)) is rounded off, and its amplitude de- 
creases with increasing strength of the scatterer potential, still vanishing with 
the system size as 1/N. The impurity lifts the degeneracy of the energy levels 
and the current then varies continuously. 8 

The model with interactions (U ^ 0) and without defect (p = 1), is solvable 
by the Bethe ansatz for periodic boundary conditions {(f) — 0) 22,23 and also for 
twisted boundary conditions (cf) ^ 0). 24-26 At half- filling, the system exhibits a 
metal-insulator transition, which occurs at U/t = 2. For U/t < 2, the system is 
in a gapless LL phase, while for U/t > 2 it is in a gapped CDW state. 27 The LL 
phase is characterized by a power-law decay of the correlations. Bosonization 
predicts that in an homogeneous LL, the leading behavior of the persistent 
current in the system size N, has a saw tooth like shape with slope —vj/(nN), 
where vj is the velocity of current excitations. 28 Since translation invariance is 
preserved in the presence of interactions, the discontinuity in the current still 
exists for finite U. 8 A Bethe ansatz calculation shows that the Drude weight 
of an homogeneous LL in the thermodynamic limit, has a finite value, which 
decreases with increasing strength of the interaction U/t. 20 The LL state is 
strongly affected by the presence of an impurity, 29-32 and bosonization yields 
that the current then vanishes as / ~ TV -1- " 3 , with a B > 0. 21 The study of 
the LL phase with p ^ 1, performed with the DMRG, 14 has in fact found this 
kind of behavior. The CDW phase is characterized by a localization length £, 
which is associated to the energy gap. From the work of Baxter, 33 the Drude 
weight in the gapped phase, is expected to behave as D <~ exp(— N/£), vanishing 
for an infinite system size. This behavior implies that although the system is 
insulating in the infinite system size limit, for a finite system, provided _/V/£ is 
not too large, D is still finite and a current can be observed. The localization 
length can then be extracted from the size dependence of the Drude weight. 

3 NUMERICAL METHOD 

We use the DMRG to numerically calculate the groundstate energy of the spin- 
less fermion system as a function of the magnetic flux, E {<p), for fixed inter- 
action U and impurity strength p, in rings up to N = 82 sites, keeping up 
to 250 density matrix eigenstates per block. 34 The DMRG is applied to the 
Hamiltonian (1) after performing the gauge transformation, which removes the 
flux into a twisted boundary condition. 11 The states of the system are char- 
acterized by the quantum numbers associated to the eigenvalues of the local 
occupation number rij and the total number of particles M = X^jli n i opera- 
tors, which commute with the Hamiltonian (l). 35 For each set of N, p and U, we 
obtained the groundstate energy E for 50 values of <fi in the periodicity interval 
— 7T < <fi < 7T, and using Chebyshev interpolation 36 we determined the corre- 
sponding current (6) and Drude weight (9), by numerical differentiation. We 
developed a DMRG algorithm for complex Hamiltonian matrices, which allowed 
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to calculate the detailed form of the persistent current / as a function of <f>, 
and to obtain the Drude weight. In a previous approach the DMRG was used 
to calculate the so called phase sensitivity AE , which is the difference of the 
groundstate energy at flux cp = and n, and can be considered a crude measure 
of the persistent current. 37-39 Although the calculation of AEo requires consid- 
erably less computational effort than the calculation of E (cf>) , because then the 
Hamiltonian matrix is real, the phase sensitivity does not provide information 
on the shape of the current and the value of the Drude weight. 

4 RESULTS 

We now present the results obtained for the persistent current and the Drude 
weight, where we take t = 1 and the interaction U is in units of t. Figs. I and 
2, exhibit NI(<fi) plotted versus <fi, for respectively, U = 0.80 and U = 3.00, 
which correspond, respectively, to the LL and the CDW phase, considering 
different impurity strengths p, on a fixed system size N = 26. We can see 
that the effect of the impurities, in both phases, is to reduce the intensity of 
the current, and to round off the shape of I(<ft). The amplitude of the current 
decreases rapidly with increasing values of \p — I|, and also with increasing 
strength of the interaction U . Fig. 3 shows that the invariance of the current 
with respect to the defect, described in Eq. (8), found for large system sizes in 
the noninteracting case, is also observed for the interacting case, both in the LL 
and in the CDW phases, the system size N required to reach the invariance being 
larger for larger interaction U . Fig. 4 displays the Drude weights associated 
to the systems with different interactions U and impurities p, fixed N = 26, 
of the currents presented in Figs. I and 2. As one would expect the Drude 
weight decreases with increasing \p — 1|, and also increasing U. Figs. 5 and 6, 
present NI{4>) plotted for several system sizes N, respectively, for U = 0.80, 
in the LL phase, and U = 3.00, in the CDW phase, with p fixed at 0.50. We 
observe that the current vanishes faster than l/N in both phases, exhibiting 
a different behavior in each phase, I{cf)) vanishing much faster with N in the 
CDW than in the LL phase. In order to analyze the behavior of the current 
in more detail we have numerically evaluated the coefficients Ik (for k = 1,2) 
of the Fourier expansion (7). The first (k — 1) and the second (k = 2) Fourier 
coefficients of the current, for U = 0.80 and U — 3.00, are shown in Fig. 
7. One can clearly see that the coefficients I\ and I2 behave similar to each 
other in both phases. However, their behavior in the LL and CDW phase is 
distinct. In the LL phase, the Fourier coefficients show a power-law decay with 
N, with the second order coefficient decaying faster, i.e., with a larger exponent, 
than the first one. In the CDW phase, the Fourier coefficients show a dominant 
exponential decay with N, with the second order coefficient also decaying faster, 
i.e., with a smaller localization length, than the first one. We observe that for 
longer rings, stronger interactions and also stronger impurities, the current is 
increasingly more precisely described by its first Fourier component. In the 
LL phase, this in fact corresponds to the asymptotic behavior predicted by 
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bosonization in the large N limit, that is / ~ Af _1_a:B sin^. However, for 
the system sizes considered here this asymptotic regime is not reached and the 
current displays a more complex behavior. The current is composed of a few 
Fourier components with decreasing weight. Fig. 8 presents the first Fourier 
coefficient of the current for different values of the interaction, U = 0.80, 2.50, 
3.00, fixed p — 0.50, from which we extract the dependence of I\ on N, in 
the intermediate range of sizes considered. We observe that for U = 0.80, in 
the LL phase, the first coefficient of the current varies as I\ ~ TV -1- " 1 , with 
ai ~ 0.06, while for U — 2.50 and U = 3.00, in the CDW phase, it varies as 
Jj ~ TV- 1 -" 5 ! exp(-iV/£i), respectively, with £i ~ 259, 6i = 0.11, and £i ~ 68, 
5\ = 0.10. The exponent ct\ is given by the slope of the straight line in Fig. 
8. a, and the length £i and the exponent 8\ were carefully adjusted in order 
to obtain the best collapse of the data in Fig.8.b, on a plot of \n(N 1+Sl I\) vs 
N/ £i . The Drude weights characterizing the systems with different interactions 
U, fixed p = 0.50, are presented in Fig. 9. Fig. 10 clearly shows that the 
results obtained for the Drude weight confirm the conducting behavior shown 
by the first coefficient of the currents in Fig. 8. We observe that, for U = 0.80 
the Drude weight varies with the system size as D ~ N~ a , with a ~ 0.04, while 
for U — 2.50 and U — 3.00 it varies as D ~ N~ s cxp(— N/£), respectively, with 
£ ~ 307, S ~ 0.08 and £ ~ 68, 5 ~ 0.06. The exponent a is given by the slope 
of the straight line in Fig. 10. a, and the localization length £ and the exponent 
S were carefully adjusted in order to obtain the best collapse of the data in 
Fig. 10. b, on a plot of ln(N°D) vs N/£. The exponents and localization lengths 
characterizing the Drude weight are a little different from those characterizing 
the first Fourier component of the current, as one would expect, since the Drude 
weight contains the contribution from the various Fourier components. One 
sees that the localization length £ and the exponent 6 decrease with increasing 
strength of the interaction U. Also, concerning the impurity influence, Fig. 4 
implies that the exponent a in the LL phase increases, and the localization 
length £ in the CDW phase decreases, with increasing \p — 1|. As mentioned, in 
the large A^ limit the current is expected to behave as its first Fourier component, 
which in the LL phase implies that the exponent a\ should be identified with 
cub = l/K — 1 as calculated from bosonization, 21 where K is the LL parameter, 
calculable from the Bethe ansatz. 40 For U = 0.80 this leads to as — 0.27, which 
is much larger than our value of a\. We should note that the size dependence 
found for the first Fourier component of the current, and the Drude weight, 
characterizes the behavior of an intermediate and limited range of system sizes. 
If one would consider a larger range of systems, in the LL phase, one would most 
probably see the data for the larger A^ bending down, crossing to an asymptotic 
power-law behaviour with the exponent approaching ag . This was observed in 
Ref. 14, where the behavior of the first few Fourier components of the current 
in the LL phase was discussed in detail, with data taken for larger values of 
A^ and stronger interaction and impurity strengths. Also, in the CDW phase 
we consider systems in an intermediate regime where the localization length 
is larger or near the system size. 41 For larger systems, the power factors that 
occur in the first Fourier component of the current and the Drude weight may 
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decline, possibly leaving a pure exponential behavior in the asymptotic regime. 

From the results obtained, we observe that the system with U = 0.80 and p = 
0.50, is characterized by an exponent a > 0, which is generated by the interplay 
of the electron interaction with the impurity, and D exhibits then a power-law 
decay with N, which implies vanishing in the limit N — > oo. On the other hand, 
the systems with U — 2.50 and U = 3.00, fixed p = 0.50, are characterized by a 
localization length £, which decreases with increasing interaction and impurity 
strength, and D exhibits now an exponential decay with N, also vanishing as 
N — > oo. Hence, we find that both in the LL and the CDW phases, with an 
impurity in the system, the effect of the interaction is to decrease the current 
and the Drude weight. As referred before, our results also provide insight into 
the spin transport in a spin-1/2 XXZ chain with a weak link, and a similar 
behavior to the one above is implied for the spin current and stiffness. So, in 
the gapless XY phase, the spin stiffness decays with a power-law, vanishing in 
the limit N — * oo, while in the gapped Ising phase, it decays exponentially, 
also vanishing as N — ► oo. Comparing our results for the persistent current 
in the LL phase, with those obtained in Ref. 14, we have similarly found that 
the current vanishes faster than 1/N. One observes that the model parameters 
strongly influence when the last asymptotic regime described by bosonization is 
reached. A calculation of the finite-size corrections to the spin stiffness in a pure 
spin-1/2 XXZ chain, 42-44 has revealed a size dependence in the gapped phase 
that has a similar form to the one found here. The result that the Drude weight 
in a ring in the gapless phase with an impurity drops to zero, is in agreement 
with a previous result obtained for a spin chain, 38 and with renormalization 
group arguments, which state that the impurity term is relevant leading to a 
transmission cut. 31,32 The renormalization group studies find that cither a weak 
barrier or a weak link lead to an insulating state for repulsive interactions, while 
in the noninteracting case those are marginal perturbations. In turn, our work 
shows that there is an invariance of the current under the defect transformation 
p — > 1/p in the interacting system, as for the noninteracting system, and that 
implies that a strong link will also reduce the current and the Drude weight. The 
observation that with an impurity in the system, the interaction always leads 
to an additional decrease of the current and the Drude weight is in agreement 
with previous results by other authors, 7 ' 8 and can be understood as it is more 
difficult to move correlated electrons in a scattering potential than independent 
electrons. 

5 SUMMARY 

We have studied the behavior of the persistent current and the Drude weight 
on a mesoscopic ring pierced by a magnetic flux. We considered a model of 
spinlcss fcrmions with repulsive interactions and a hopping impurity, which is 
also equivalent to a spin-1/2 XXZ chain with a weakened exchange coupling. 
Using a powerful numerical method, the DMRG with complex fields, we have 
calculated the detailed form of the current as a function of the magnetic flux, 
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which enabled us to investigate the corrections to the large system-size limit, 
and also allowed to obtain the Drude weight. We show that the system at half- 
filling, changes from an algebraic to an exponential behavior as the interaction 
increases, corresponding to a change from a LL to a CDW phase. We find that 
the analytical predictions of bosonization for the LL phase, are not yet fully 
observed in the intermediate range of system sizes considered. In addition we 
observe that the invariance of the current under the defect transformation p — > 
1/p, seen in the noninteracting system, is also verified in the interacting system, 
in both phases. Hence, an isolated strong link is not only useless for increasing 
the persistent current (as might have been expected), but it rather destroys 
coherence and reduces the current. The behavior determined for the current is 
consistent with the behaviour determined for the Drude weight, the LL phase 
being characterized by an exponent a > 0, which results from the interplay of 
the interactions with the impurity, while the CDW phase is characterized by a 
localization length £, which decreases with increasing interaction and impurity 
strength. We find that, both in the LL and the CDW phase, with a defect in the 
system the interactions always suppress the current, and the Drude weight drops 
to zero in the limit N — > oo. Away from half-filling there is no metal-insulator 
transition in the pure case, and the system is always metallic. Hence, one 
does not expect to observe then a change in the current and the Drude weight 
from an algebraic to an exponential decay. Nevertheless, one still expects to 
observe that in the system with an impurity, the current and the Drude weight 
decrease with increasing impurity and interaction strcngthes. Therefore, within 
the model considered, the interactions cannot explain the results observed in 
the experiments. 



*Electronic address: fdias@cii.fc.ul.pt 

References: 

1 L.P. Levy, G. Dolan, J. Dunsmuir, and H. Bouchiat, Phys. Rev. Lett. 64, 
2074 (1990). 

2 V. Chandrasckhar, R.A. Webb, M.J. Brady, M.B. Ketchen, W.J. Gallagher, 
and A. Kleinsasser, Phys. Rev. Lett. 67, 3578 (1991). 

3 D. Mailly, C. Chapelier, and A. Bcnoit, Phys. Rev. Lett. 70, 2020 (1993). 
4 M. Biittiker, Y. Imry, and R. Landauer, Phys. Lett. 96A, 365 (1983). 
5 Y. Imry, Introduction to Mesoscopic Physics, Oxford University Press (1997). 
6 U. Eckern and P. Schwab, Adv. Phys. 44, 387 (1985); J. of Low Temp. Physics 
126, 1291 (2002). 

7 M. Abraham and R. Bcrkovits, Phys. Rev. Lett. 70, 1509 (1983). 

8 G. Bouzerar, D. Poilblanc, and G. Montambaux, Phys. Rev. B 49, 8258 

(1994). 

9 S.R. White, Phys. Rev. Lett. 69, 2863 (1992). 

10 I. Peschel, X. Wang, M. Kaulke, and K. Hallberg (Eds.), Density- Matrix 
Renormalization - A New Numerical Method in Physics, Springer (1999). 
n U. Schollwock, Rev. Mod. Phys. 77, 259 (2005). 



9 



12 J. Voit, Rep. Prog. Phys. 58, 977 (1995). 

13 See, e.g., G. Griiner, Density Waves in Solids, Addison- Wesley Publishing 
Company (1994). 

14 V. Meden and U. Schollwock, Phys. Rev. B 67, 035106 (2003). 

15 See, e.g., G. D. Mahan, Many-Particle Physics, Plenum Press (1990). 

16 N. Byers and C.N.Yang, Phys. Rev. Lett. 7, 46 (1961). 

17 M. Henkel and D. Karcvski, Eur. Phys. J. B 5, 787 (1998). 

18 M. Henkel, Conformal invariance and Critical Phenomena, Springer (1999). 

19 W. Kohn, Phys. Rev. 133, A171 (1964). 

20 B.S. Shastry and B. Sutherland, Phys. Rev. Lett. 65, 243 (1990). 

21 A.O. Gogolin and N.V. Prokof'ev, Phys. Rev. B 50, 4921 (1994). 

22 C.N. Yang and CP. Yang, Phys. Rev. 150, 321 (1966). 

23 J. Des Cloizeaux and Michel Gaudin, J. Math. Phys. 7, 1384 (1966). 

24 C.J. Hamer, G.R.W. Quispcl, and M.T. Batchclor, J. Phys. A 20, 5677 (1987). 

25 F.C. Alcaraz, M.N. Barber, and M.T. Batchclor, Phys. Rev. Lett. 58, 771 

(1987); Ann. Phys. 182, 280 (1988). 

26 B. Sutherland and B.S. Shastry, Phys. Rev. Lett. 65, 1833 (1990). 

27 In the XXZ chain, a transition occurrs at A/J = 1, from an XY model, for 

A/J < 1, to an Ising model, for A/J > 1, A/J = 1 corresponding to the 

Hcisenberg model. A gap opens up in the spin-excitation spectrum, that is 

equivalent to the gap in the charge-excitation spectrum of the spinless fermion 

system. 

28 F.D.M. Haldanc, J.Phys. C 14, 2585 (1981). 

29 W. Apel and T.M. Rice, Phys. Rev. B 26, R7063 (1982). 

30 T. Giamarchi and H.J. Schulz, Phys. Rev. B 37, 325 (1988). 

31 C.L.Kanc and M.P.A. Fisher, Phys. Rev. B 46, 15233 (1992). 

32 S.Eggert and I. Affleck, Phys. Rev. B 46, 10866 (1992). 

33 R.J. Baxter, Exactly Solved Models in Statistical Mechanics, Academic Press, 

(1982). 

34 The number m of eigenstates required to achieve convergence, of 1 part in 10 5 , 
for the ground state energy E , increases with the system size, e.g., m = 128 
for N = 30 and m = 250 for N = 66. 

35 For the spin XXZ chain the quantum numbers characterizing the states cor- 
respond to the local spin projection S? and the total spin projection S z = 

36 W.H. Press, S.A. Teukolsky, W.T. Vetterling, and B.P. Flannery, Numerical 
Recipies in Fortran 77 - The Art of Scientific Computing (Second Edition), 
Cambridge University Press (2001). 

37 P. Schmitteckert, T. Schulze, C. Schuster, P. Schwab, and U. Eckern, Phys. 
Rev. Lett. 80, 560 (1998). 

38 D. Weinmann, P. Schmitteckert, R.A. Jalabert, and J.L. Pichard, Eur. Phys. 
J. B 19, 139 (2001). 

39 T.M.R. Byrnes, R.J. Bursill, H.-P. Eckle, C.J. Hamer and A.W. Sandvik, 

Phys. Rev. B 66, 195313 (2002). 

40 F.D.M. Haldane, Phys. Rev. Lett. 45, 1358 (1980). 



10 



41 N. Laflorencie and H. Rieger, Eur. Phys. J. B 40, 201 (2004). 
42 S-J. Gu, V.M. Pcrcira and N.M.R. Peres, Phys. Rev. B 66, 235108-1, (2002). 
43 N. Laflorencie, S. Caponni, and E.S. S0rensen, Eur. Phys. J. B 24, 77 (2001). 
44 F. Heidrich-Meisner, A. Honecker, D.C. Cabra, and W. Brenig, Phys. Rev. B 
68, 134436 (2003). 



11 



Figure Captions: 

FIG. 1. Persistent current NI(<j>) vs 0, at fixed N = 26, for U = 0.80 and 
different p. 

FIG. 2. The same as in Fig. 1, but for U = 3.00. 

FIG. 3. Symmetry of the current NI(<p) in p (empty symbols) vs 1/p (filled 
symbols), for U = 0.80 (circles) at N — 18, and U = 3.00 (diamonds) at 
N = 58. 

FIG. 4. The Drude weigth D vs p, at fixed TV = 26, for U = 0.80 and U = 3.00. 
FIG. 5. Persistent current NI(<f>) vs <j>, for U = 0.80, p = 0.50 and increasing 
N. 

FIG. 6. The same as in Fig. 5, but for U = 3.00. 

FIG. 7. Fourier coefficients of the current. \n(NIk) vs ln(N), for U = 0.80 
(circles) and U = 3.00 (diamonds), at p = 0.50, for k = 1 (filled symbols) and 
fc = 2 (empty symbols). 

FIG. 8. Nh for different U and p = 0.50. (a) ln(JV7i) vs ln(N), for f7 = 0.80, 
the line represents h - iV" 1 -" 1 . (b) ln(iVJi) vs N, for U = 2.50 and 3.00, the 
lines represent I\ ~ iV _1_<51 exp(— N/£i), with and (5i dependent on U. 
FIG. 9. The Drude weigth Uasa function of 1/N, for different U and p = 0.50. 
FIG. 10. D for different U and p = 0.50. (a) ln(L>) vs ln(N), for [/ = 0.80, 
the line represents D - A^ Q ; (b) ln(/J) vs iV, for U = 2.50 and 3.00, the lines 
represent D ~ iV _<5 exp(— iV/^), with ^ and 5 dependent on U. 
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